Population dynamics in a random environment 
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We investigate the competition between barrier slowing down and proliferation induced superdif- 
fusion in a model of population dynamics in a random force field. Numerical results in d = 1 suggest 
that a new intermediate diffusion behaviour appears. We introduce the idea of proliferation assisted 
barrier crossing and give a Flory like argument to understand qualitatively this non trivial diffusive 
behaviour. A one-loop RG analysis close to the critical dimension dc = 2 confirms that the random 
force fixed point is unstable and flows towards an uncontrolled strong coupling regime. 



(N 



> 
oo 

o 
o 
o 



o 
o 



The presence of disorder often radically changes the 
statistical properties of random walks. For example, ran- 
dom walks in a random potential are trapped in deep 
potential wells: this may lead to sub- diffusion, i.e. the 
fact that the typical distance traveled by the walkers 
grows more slowly than the square-root of time A 
much studied model exhibiting this type of behaviour 
is the Sinai model, where particles diffuse in a random 
force field in one dimension In this case, the en- 

ergy barriers typically grow as the square root of the 
distance, which leads to a logarithmically slow progres- 
sion of the random walkers. There are also several mech- 
anisms which lead to .super- diffusion. For example, if 
the random force field is rotational, the random walkers 
can be convected far away by long streamlines An- 
other interesting mechanism of superdiffusion is random 
proliferation: suppose that each random walker can ei- 
ther die or give birth to new random walkers at a rate 
which is random, both in time and space. There is in 
this case a possibility for an 'outlier' random walker, 
that has by chance traveled a distance much greater than 
the square-root of time, to have been particularly pro- 
lific: he and his siblings then represent an appreciable 
fraction of the whole population, leading to a motion 
of the center of mass faster than diffusive. This mech- 
anism has been much studied (although not explicitely 
discussed as such) in the context of Directed Polymers 
(dp) in random media or equivalently the Kardar-Parisi- 
Zhang (kpz) model of surface growth ||]. The aim of 
the present work is to investigate the case where both 
these mechanisms are present simultaneously. The mo- 
tivations for such a mixed model are numerous. In the 
context of population dynamics (for example, bacteria on 
a random substrate), similar models have recently been 
investigated, with quite interesting results One can 
also give an economic interpretation of population dy- 
namics, where the local density of random walkers is the 
wealth of a given individual. Biased diffusion represents 



trading between individuals, whereas the random growth 
term is the result of speculation |^ . One can argue that 
generically, this type of model leads to a Pareto (power- 
law) tail in the distribution of wealth . Finally, from a 
theoretical point of view, this mixed model leads to the 
interesting possibility of new behaviour, intermediate be- 
tween superdiffusion and subdiffusion. 

More precisely, we study here the following equation 
for the local population density P{x, t) in d dimensions: 



dP{x, t) 
di 



j/oAP(x, t) - V{F{x)P) + ?7(f , t)P{x, t), 



(1) 



where vq is the bare diffusion constant, F(x) a space 
dependent static Gaussian random force such that 
{F^{x)F,{x'))f = 2cj15^,.,5\x - x') §, and 7?(a-,t) a 
Gaussian random growth rate, depending both on space 
and time, with {r]{x,t)r){x' ,t'))^ = 2al5{t-t')5'^{x-x'). 
Due to the last term, the total population Z{t) — 
J dxP{x,t) is not conserved. The quantities of interest, 
which describe how the population spreads in time are, 
for example, the average center of mass motion, 



— J xP{x, t) dx 



/F,ri 



(2) 



or the average 'width' of the diffusing packet (A^) : 



xP{x, t) dx (3) 



A2(i) = ^ / x'^P{x,t) dx 



(Other moments can however also be studied: see below). 
An alternative description is in terms of the free-energy 
h(x,t) — log P{x,t), which obeys the equation: 



dh{x, t) 
dt 



uoAhix, t) + X{Vhf 

F{x) • V/i - V • F{x) + ri{x, t), (4) 
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with \ = i/Q. When F = 0, these equations represent the 
well-known kpz (or Directed Polymer) problem, whereas 
for 77 = 0, one recovers the problem of a random walk 
in a random environment. Both problems can be ap- 
proached using a perturbative Renormalization group; 
interestingly, the critical dimension for both problems is 
dc — 2. For the random drift problem, one finds that 
the coupling constant gp = ap / {2TT)h'Q flows towards a 
non trivial fixed point of order e in dimensions d = 2 — e 
P,p^P]. This in turn leads to a subdiffusive behaviour: 
Xcmit) grows as tp with j/^ = (1 — e^)/2 < 1/2. 
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FIG. 1. One loop RG flow in the gn,gF plane. As soon as 
Qn is non zero, it flows towards the strong coupling region. 
Our numerical simulations suggest that an attractive 'mixed' 
fixed point appears in this region. 

For the kpz problem, the coupling constant is = 
a"^}? / {2-k)vq] the Gaussian fixed point = is again 
unstable for c? < 2, but there is no accessible fixed point 
at one loop for d > 3/2 The exponent v therefore 

cannot be computed but is expected (and found numeri- 
cally) to be greater than 1/2: in the population dynamics 
language, the possibility of far-away proliferation leads to 
superdiffusion. On the other hand, for d < 3/2, a physi- 
cal fixed point appears; in d = 1, the one loop calculation 
even provides the exact result Vri = 2/3 for reasons de- 
tailed in |11|. We have performed a RG analysis in the 
mixed case where both gp and g^, are non zero. This can 
be done using a field theoretical representation (Martin- 
Siggia-Rose) representation of Eq. (^ , which allows one 
to generate the perturbation expansion in gp and g,,. 
Performing calculations along the lines of ]9|-pl|] , we find 
that the two (3 functions are given by: 



dgr, 

dl 



e.9r, + 2.g,,.g_F 
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dgp_ 

dl 



9f, 



(5) 



where £ is the logarithm of the running length scale. The 
resulting fiow is represented in Fig. 1. The subdiffusive 
random force fixed point g,i = 0, gp — e is therefore un- 
stable in the presence of a small 'proliferation' term 5^. 
Unfortunately, at one loop, fiows towards the strong 
coupling region, as is the case in the standard kpz case 



gp = 0. 

In order to obtain some information about this strong 
coupling behaviour, we have performed some numerical 
simulations in one dimension, where both the fixed points 
corresponding to Sinai subdiffusion and to DP/kpz su- 
perdiffusion are rather well understood. We have found 
results that suggest the existence of an attractive fixed 
point, characterized by a new non trivial diffusive be- 
haviour (intermediate between the Sinai and dp / kpz be- 
haviour) . We have numerically evolved a space and time 
discretized version of Eq. (|^), and have worked with 
logP to avoid precision problems. Starting from a lo- 
calized packet P{x = ia,t — 0) — Si,Q, we have found 
that as soon as both coupling constants gp and 5° are 
non zero, the exponent v describing the diffusion of the 
center of mass Xcm{t) at large times is found to be close 
to the value v* = 1/2 (see Fig. 2). The position Xmax{t) 
of the maximum of P{x, t) behaves very similarly. 
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FIG. 2. Behaviour of the average center of m&SS "Xcm (t) and 
of the average position of the maximum of the packet Xmax (t) 
(i.e. the point where P{x,t) is maximum), as a function of 
time, for a-qjap = 0.125. The best linear fits are shown, and 
lead to the an estimate for slightly smaller than 1/2. Inset: 
Value of 1/* (determined from the behaviour of Xcmit)) as a 
function of the Hurst exponent of the potential H, compared 
with the Flory prediction. 

The ratio gp/g^ affects only the short time transient 
behaviour, which is either Sinai like or dp/kpz like, as 
shown in Fig. 3. In the RG language, this suggests that 
a non trivial attractive fixed point gp,g*j appears. This 
is compatible with the fiow diagram of Fig. 1 ||l^, al- 
though this new fixed point is out of reach at the one- 
loop level. Although the value i/* = 1/2 corresponds 
to free-diffusion, the motion of the packet for a given en- 
vironment is very far from simple diffusion, as the study 
of the width A of the packet shows. We have found nu- 
merically that {A'^)p^r] behaves as t'^''" with Ci — 0.24, 



2 



C2 ~ 0.34 and (4 ~ 0.38. This non trivial behaviour 
is actually present for both the Sinai problem and the 
DP/kpz problem. This comes from the fact that for both 
problems, the effective free energy h{x,t) = log P{x,t) 
behaves as a random walk in x space. This is triv- 
ial for the Sinai problem, since the potential is indeed 
constructed as the sum of local random forces. For the 
dp/kpz problem, this is far less trivial and results from 
the fact that one can obtain exactly the stationary dis- 
tribution of h{x,t) in one dimension, which turns out to 
be the same as for the linear case A = 0, i.e., again a 
random walk in x space It is well known that for a 
random walk potential, the probability that two nearly 
degenerate minima are separated by a distance A falls 
off as A^^/^ for large A. The qth moment of A is there- 
fore dominated by extreme events as soon as q > 1/2. 
Physically, this means that for most realizations of the 
disorder, the width A of the packet is small ||l3| , |r^ , ex- 
cept in rare situations where the packet is divided into 
two subpackets very distant from one another. The natu- 
ral cut-off for A is of the order of Xcm {t) itself. Therefore 
one obtains, for q > 1/2, (A'?(t)) cx [xcmit)]''''^/'^ ■ For 
the Sinai problem, using Xcmit) oc \og^{t), this leads to 
(A^(t))i? oc log^{t), whereas for the dp/kpz case, using 
Xcm{t) oc t^/^, one finds (A^(t))^ cx t: both these results 
are actually exact, as has been shown in fi p^JTs] . As- 
suming that the effective potential in the mixed case is 
again a random walk in x space, and using v* ~ 1/2, we 
obtain C, = (2q - l)/4(j, i.e. Ci ^ 0.25, C2 ^ 0.375 and 
(■4 ~ 0.4375, in reasonable agreement with our numerical 
values p^ ]. 
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FIG. 3. Behaviour of the average center of mass Xcm,{t) for 
different values of ap- This curves shows that for large ctf, 
the short time behaviour is Sinai like, crossing over to the 
mixed behaviour at large times. The value of Xc,n{t) has been 
rescaled by a^'^ to obtain a reasonable data collapse at large 
times. 



In order to understand the value of z^* ~ 1/2, one 
needs to develop a consistent picture of the competition 
between the slowing down induced by the ever-growing 
Sinai barriers and the speeding up of the population 
spreading allowed by the multiplicative growth term 77. 
Before addressing the full Sinai-|-KPZ problem, we first 
consider the simpler case of a unique barrier of height 
Uq, which develops on scale L. For definiteness, we have 
solved numerically the equation (|I]) on the interval [0, L] 
with F{x) = —Uo/Ls'm{ATrx/L). The initial condition is 
localized in the first well, and the crossing time r is de- 
fined as the average time after which the relative weight 
of the population in the second well is half of that in the 
first well. For rj = 0, one finds the classical Arrhenius 
law: logr = Uq/i^o- When 77 7^ 0, the behaviour of r as 
a function of Uq for different values of L is shown in Fig. 
4. The result can be expressed as: t oc L"^^^ f{Uo/ VL), 
with f{y 0) = 1 and f{y — > 00) oc y''. This scaling of 
r with L can easily be understood. In the limit Uo — *■ 0, 
the time for the particles to reach a distance L is given 
by the dp/kpz scaling, i.e. L oc r^^^. 
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FIG. 4. Average barrier crossing time (r), rescaled by L^^^, 
as a function of the barrier height Uo rescaled by y/L, for dif- 
ferent sizes L, and in log-log coordinates. The slope 5/3 is 
shown for comparison. The power law increase of (r) as a 
function of Uo has to be compared to the usual exponential 
(activated) dependence. 

The influence of the external potential Uq becomes 
substantial when it becomes of the order of the effec- 
tive KPZ potential h, which, as discussed above, grows 
as Numerically, the exponent b is found to be very 
close to 6 = 5/3. Therefore, the exponential increase 
of the crossing time with the barrier height is replaced 
by a power-law increase in the presence of the random 
growth term rj. One can call this effect proliferation 
assisted barrier crossing: the probability that a parti- 
cle reaches the top of the barrier x* by pure diffusion 



3 



is exp(— f/o/i'o); but due to the random growth term, 
this probabihty is multipHed by a certain prohferation 
'gain' factor ex.pQ{x* ,t) If the path C leading from 
the initial point xq to x* was unique, one would simply 
have G{x* , t) — J dt'ri{xc{t'), t'), which typically behaves 
as (T,j\/t. In fact, many paths contribute to Q{x*,t). 
This leads to a kind of preaveraging effect of the ran- 
dom growth term rj over the width w{t) of the paths C. 
Therefore: 



dt' 



w{t 



1/2 



(6) 



Since most paths leading to x* spend their time in the 
thermally accessible region of the well, one can estimate 
w{t') as w = L/\/Uo. The proliferation factor then com- 
pensates the barrier when t oc Uq^"^ (for d = 1). This 
simple argument therefore leads to b = 3/2, not very far 
from the numerical value b ~ 5/3. Actually, one can ap- 
ply this argument to the unconfined case Uq = 0, where 
the detrimental factor is now the entropy of the random 
walk exp{—x*^/t). Using self-consistently 'w{t') — x*{t'), 
the compensation argument now leads to x* cx i3/(4+(i)^ 
which is precisely the Flory result for the DP / KPZ prob- 
lem. This Flory value can be obtained using a variational 
method, either with replicas or without replicas [^ . 
In spirit, Eq. (||) is actually very close to the latter cal- 
culation. The value & = 3/2 can therefore be seen as a 
Flory value for this problem. 

Returning now to the Sinai case, where the barrier 
height grows as ap^/x* , the self consistent compensa- 
tion argument now leads to ap^/x* ~ ari\/t/x* ^ or 
X* ^ (unl OF)\Jt- The \ft behaviour is close to the nu- 
merical results shown in Fig. 2. However, as shown in 
Fig. 3, the dependence of Xcm on ct_f is found to be weaker 
than the I/ctf behaviour predicted by this simple argu- 
ment, and closer to l/cr^''*. We have also investigated 
numerically the case where the force derives from a frac- 
tional Brownian motion with a Hurst exponent H . The 
case H ~ 1/2 is the standard Sinai random walk poten- 
tial considered above. An extension of the proliferation 
argument to this case predicts that v* = 1/(1 -f 2iJ) for 
H > 1/4, reverting to the dp/kpz value v* — 2/3 for 
smaller values of H (i.e., when the potential is not 'con- 
fining' enough). As shown in Fig. 2, our numerical values 
for ly* agree quite well with this prediction: for example 
p*{H = 3/4) ~ 0.37 and v*{H = 1/4) ~ 0.65. 

In summary, we have investigated the competition be- 
tween barrier slowing down and proliferation induced su- 
perdiffusion in a model of population dynamics in a ran- 
dom force field. The one-loop RG analysis close to the 
critical dimension dc ~ 2 predicts that the subdiffusive 
fixed point is unstable against 'proliferation' and flows 



to strong coupling. Our numerical results in d = 1 ac- 
tually suggest that both the Sinai and kpz fixed points 
are unstable, and flow towards a new stable mixed fixed 
point. We have given a heuristic Flory like argument, 
which allows us to understand qualitatively the diffusive 
behaviour at this mixed fixed point, and also our results 
on proliferation assisted barrier crossing. This work can 
be extended in various directions: for example, a two- 
loop RG calculation would be interesting. One could also 
study the effect of non linear terms in the population 
equation, such as —P^ or V • (PVP), and the role of a 
non zero external force {F{x)). It would be worth per- 
forming some numerical simulations of the barrier cross- 
ing problem and of the mixed model in c? = 2. 
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